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Abstract 

The distribution of individual Dirac eigenvalues is derived by relating them to the density and 
higher eigenvalue correlation functions. The relations are general and hold for any gauge theory 
coupled to fermions under certain conditions which are stated. As a special case, we give examples 
of the lowest-lying eigenvalue distributions for QCD-like gauge theories without making use of 
earlier results based on the relation to Random Matrix Theory. 
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1 Introduction 



Consider the partition function of Yang-Mills theory coupled vector-like to fermionic matter in a given 
representation of the gauge group. Integrating out the Nf fermions in the Euclidean partition function 
one obtains the formal expression 



Z = J[dA] Y[det[0(A)-m f ] exp[-S[A]] , (1.1) 



where S[A] is the pure gauge theory action. This involves the determinant of the Dirac operator L/> 
(plus the mass terms). One way to make the evaluation of this determinant meaningful is to regularize 
the theory in such a way that there is a finite number of Dirac operator eigenvalues Aj. 

A regularization of the theory in such a way already exists: lattice gauge theory with a finite 
lattice spacing a and a finite four-volume V. Because of chiral symmetry and the index theorem the 
spectrum of the Dirac operator eigenvalues is constrained: the spectral density p(A) will be symmetric 
around the origin, p(X) = p(—X), and there will, for any given gauge field configuration, be a number 
v of exact zero modes. This number of zero modes is dictated by the index theorem. 

What will be the distribution of the remaining eigenvalues? The spectral density gives the overall 
distribution of the unordered set of eigenvalues, but one would like to ask more detailed questions. 
Here the lattice gauge theory approach offers a fruitful angle to the problem. It is thus instructive to 
view the path integral over the gauge field A^x) as the limit of a discrete sum over gauge field config- 
urations. For each gauge field configuration one will have a set of non-zero Dirac operator eigenvalues 
Aj, where, because of the spectral symmetry, only the N positive ones need to be considered here. 
These eigenvalues can be ordered according to their magnitudes, Ai < A2 < • • • < Xn- Performing 
this ordering configuration by configuration, one can ask for the individual distribution of the kth 
Dirac operator eigenvalue. It is far from obvious how one would go about computing such individual 
distributions of ordered Dirac operator eigenvalues directly from first principles and the quantized 
gauge field theory. In fact, it is surprising how little apparently can be said from first principles. The 
only rigorous statement we have been able to find in the literature is the inequality 

Ai < j , (1-2) 

where Ai is the smallest non-zero eigenvalue, C is a constant, and L is the linear extent of the four- 
volume. All this says is that the smallest eigenvalue Ai approaches zero in the limit L — > 00 at a rate 
which is not slower than that of free fields. This should be contrasted with what one obtains if one in 
addition assumes spontaneous chiral symmetry breaking. From the Banks-Casher relation (2j 

p(0) = L 4 £/vr , (1.3) 

where £ is the chiral condensate, it follows that a non-vanishing £ requires the Dirac operator eigen- 
values to accumulate at substantially faster rate, yielding (with some other constant C) |Sj 

Ai < § • (1-4) 

This has been succinctly rephrased by Shuryak and Verbaarschot [3]: Spontaneous chiral symmetry 
breaking and the required accumulation of eigenvalues near the origin according to (|1.4j) suggests that 
a more appropriate object to study near A ~ is the so-called microscopic spectral density 

P.(0 = ±p(x = £ 7 ) . (1-5) 
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Remarkably, it has turned out that this microscopic spectral density is computable analytically using 
the relation to the appropriate low-energy effective field theory, a chiral Lagrangian What is per- 
haps more astonishing, the same results can be obtained from universality classes of large- N (chiral) 
Random Matrix Theories [U 113 El IE]- To be precise, using the relation to effective field theory it has 
been shown how to derive the microscopic spectral density p s (C) M an d spectral two-point function [§| 
for the symmetry breaking class of QCD. For the two remaining symmetry breaking classes the corre- 
sponding effective field theories have been compared perturbatively to the microscopic spectral density 
from chiral Random Matrix Theory ^HJ- For the symmetry breaking class of QCD the procedure for 
deriving the higher /c-point spectral correlation functions using effective field theory is completely 
clear, and it is only computational tedium that has so far prevented the calculation of higher fc-point 
functions. Very compact expressions for these spectral /c-point functions are known if one accepts the 
conjecture of the relation to Random Matrix Theory, and the relevant generating functional which 
eventually should be derived from field theory is also known explicitly jlll I12j . But in this paper we 
will not rely on chiral Random Matrix Theory at all, and therefore we shall not make use of these 
explicit results, even though they are almost certainly in agreement in what one will obtain from 
effective field theory. 

The purpose of this paper is to point out that even the distribution of individual Dirac operator 
eigenvalues can be computed directly from effective field theory. This may seem surprising, because 
the concept of an ordered set of Dirac operator eigenvalues is not obviously related to quantum field 
theory observables. But following the program described in refs. , it is clear how to proceed if one 
instead wishes to derive all spectral correlation functions. Having all spectral correlation functions 
there should be no more spectral information in the theory, and individual eigenvalue distributions 
ought to follow. This is indeed the case. In fact, the derivation will be completely general, and 
can be phrased in terms of fundamental notions from probability theory and certain assumptions 
regarding the regularized QCD partition function which we will state explicitly below. In particular, 
the discussion will not be restricted to the "e-regime" of QCD, although this clearly is an important 
special case. The counting of eigenvalues from the origin can trivially be replaced by the counting 
from any other arbitrary point in the spectrum. 

2 Spectral Correlation Functions and Eigenvalue Distributions 

To proceed, we need precise definitions of the spectral correlation functions. We remind the reader 
that we consider a situation in which we have a finite number of Dirac operator eigenvalues, and we 
focus on the N eigenvalues lying above zero. To even define spectral correlation functions of these iV 
Dirac operator eigenvalues Aj=i jv we must assume the existence of a joint probability distribution 
function Vn(\i, ■ ■ ■ , Ajv). For simplicity we will here assume that this function has support on the 
whole space of (the restriction to a subset hereof is immaterial for the discussion that follows, and 
is readily implemented). We will also assume that Pjv(Ai, . . . , Ajy) is symmetric under interchange of 
all arguments. We choose the normalization convention 



The limit N — > oo is easily taken in all formulas that follow below. 

The /c-point density correlation function, which gives the probability density to find the k variables 
at values Ai 3 . . . , A& is then 




(2.1) 




(2.2) 
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In quantum field theory one is most used to considering just the one-point function, the spectral density 
p(A). As follows, it is obtained by integrating the joint probability density over N — 1 variables, viz. 



p(A) = N d\ 2 ...d\ N V N (\,\ 2 ,...,\ N ) , (2.3) 
J o 

the factor of N easily being understood as due to having arbitrarily chosen to keep just the first entry 
fixed while integrating out the rest. The existence of a joint probability density Vn{A\, ■ ■ ■ , An) is 
implicit in all work that considers the spectral density and higher spectral correlation functions of 
the Dirac operator in QCD and related field theories. In particular, the iV-point function is simply 
proportional to the joint probability distribution function itself, pn(Ai, . . . , An) = Nl Vn(Ai, . . . , An)- 
It is useful to also define a fc-th "gap probability". For simplicity, let us restrict ourselves to a 
gap adjacent to the origin; it is straightforward to pick any other fixed interval on the support of the 
eigenvalues. Define 



Ek{3) 55 JN=Q 



ps poo 

-J dA x ...dA k J dA k+1 ...dA N T N (A 1 ,...,A N ) , for k = 0, 1, . . . ,N . (2.4) 



For k = this gives the probability that all variables are larger or equal to s, Aj = i ... jv > s, meaning 
that the interval [0, s] is free of eigenvalues. For general k the interval [0, s] is occupied by k eigenvalues 
and [s, oo) is occupied by N — k eigenvalues. Our first step consists in proving that the k-th gap 
probability can be written as a sum of integrals over (N — k + 1) different /-point functions. We use 
the simple identity 

i / -\ 



(a -by = ^(-r/M U-^ , (2.5) 



and choose a = dA and b = Jq dA to replace all the (N — k) integrals dA in eq. (|2.4jl by a — b: 



■V 



TV' 

E k (s) = dAx.-.dA, 



fN — k\ ( f'°°\ JV_ft_( / f s \ 1 
^ (_1) V 1 [JJ dA k+1 ...dA N V N (A 1 ,...,AN) 

N-k ^ „ s 

(-I)'tj / dAt . . . dA k+l Pk+l (Ai,. . . , A k+l ) . (2.6) 
i=o ' ^° 

Here we have used the invariance of Vn(A\, ■ ■ ■ , An) under permutations. The formula (|2.6|) neatly 
expresses the gap probability in terms of spectral correlation functions For the simplest case of 
k = we get 

E (s) = 1 - l dAipi(Ax) + ~ / dAidA 2/92 (Ai,A 2 ) -\ f dA 1 dA 2 dA 3 p 3 (A 1 ,A 2 ,A3)+ ••• • (2-7) 
Jo 2 Jo 6 Jo 

We may also introduce a generating functional for all gap probabilities, defined as |13j 

N 1 ^ 

E(s; 6 = 1 + E(-^7i / dXl ■ ■ ■ dXl • • • , A ; ) . (2.8) 

i=i ' Jo 



It immediately follows that 



, for k = 0,l,...,N . (2.9) 

£=1 
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The probability to find the k-th eigenvalue at value Xk = s is 



N 



p k (s) = k [^ k J J d\i . . . dXk-i J dX k +i ■ ■ .d\ N T>N{Xi, . . . , X k -i, Xk = s, Xk+i, ■ ■ ■ , Ajv) ■ (2.10) 

If we order the eigenvalues, Ai < . . . < AjV; the quantity Pk(s) gives the probability to find the fc-th 
eigenvalue in this ordering at Xk = s, as ea. H2.1U|) gives the probability to find k — 1 eigenvalues in 
[0, s], Xk = s, and N — k eigenvalues in the complement [s,oo). It is shown in Appendix 1X1 that 
the distributions of eigenvalues Pk(s) are indeed all properly normalized, J °° ds Pk(s) = 1 for all 
k = l,...,N. 

If we compare eq. (|2.1()|) and eq. (|2,4j) for k — 1 we observe that they differ only by one integration. 
By differentiating we thus easily obtain 

^E k {s) = k\ (p k {s) - Pk+1 (s)) , (2.11) 

defining po(s) = 0. Here we have again made use of the permutation symmetry of the integrand, 
renaming the variable of the differentiated integrals by A^. Inserting eq. (|2.6f) or eq. ()2.9|) we can 
thus recursively determine the distribution of the k-th individual eigenvalue from the knowledge of all 
density correlation functions alone. This confirms the intuition mentioned in the introduction that if 
all spectral correlations are known, no more spectral information can be gleaned. 
Let us give some examples. For the first eigenvalue we obtain 

Pi(s) = ~§- s E o(s) = Pi(s)- J* dXp 2 (X,s) +... , (2.12) 

by a single differentiation of eq. (|2.7|) . whereas the second eigenvalue is given by 

p 2 {s) = [ dX P2 {X,s) + ... . (2.13) 



J o 

This explicitly confirms the expectation that the first approximation to the smallest of these ordered 
eigenvalues is given by the spectral density itself. The higher order correlation functions systematically 
correct for the error in this initial approximation. As we shall show below for the case of the microscopic 
spectrum of the Dirac operator, the convergence can be remarkably fast. Proceeding recursively, an 
infinite sequence of eigenvalue distributions follows from eq. Q2.11[) . 



3 Distributions of Dirac Operator Eigenvalues in the e-Regime 

An obvious place to apply the above general formalism is near the origin. As mentioned in the 
introduction, if one assumes that chiral symmetry is broken spontaneously the spectral correlation 
functions for the lowest-lying Dirac operator eigenvalues can be computed analytically. The bridge 
is an effective field theory which hinges on the precise Goldstone boson manifold, and the pertinent 
field theory is taken to a corner known as the e-regime |14j by a suitable tuning of the four-volume 
V. Let us mention that if we want to compare the regularized gauge theory with a finite number of 
eigenvalues N to such an effective theory we must assume that this number is large. 

The density correlation functions and thus also the distributions of eigenvalues depend on the 
number of flavors Nf, quark masses rrif and gauge field topology v. There have been several numerical 
simulations of individual Dirac eigenvalue distributions both using staggered fermions (see, e.g., refs. 
|151 I16j ) and overlap fermions [171 1181 I19j . So far it is has only been possible to compare with 
the analytical predictions based on the conjectured Random Matrix Theory results |20| I21j. and an 
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outstanding question has been whether these individual eigenvalue distributions also can be derived 
directly from the effective field theory. We will here rely only on what to date has been derived from 
field theory, namely the one- and two-point functions for the symmetry breaking class of QCD-like 
gauge theories (gauge theories with fermions transforming according to complex representations of the 
gauge group) [3J |§] . 




Figure 1: The distribution of the first Dirac eigenvalue for Nf = and \i>\ = 0,1,2. The first term 
(green) and the sum of first and second term (red) in the expansion eq. (|2.12|) are compared to the 
full result from Random Matrix Theory (blue) summing over all density correlations. 

For simplicity of illustration we restrict ourselves to massless quarks. The microscopic spectral 
density eq. (|1.5|) for Nf massless flavors in a sector of topological charge v derived in [3] reads 

pl N " v \0 = ^[4 /+H (C)-Jjv /+ M + i(C)JiV /+ H-i(C)] , (3-1) 

in terms of the rescaled eigenvalues £ = XVT,. The two-point function has been derived so far only 
for zero flavors Nf = and arbitrary v ^ 

P { s Nf=0 ' u \Cv) = P f' u \0pf^(r,) - ^ [CJ M+1 (C)J M (iy) - H^) J H(0] 2 • (3-2) 

We can now give the approximate distribution of the first Dirac Eigenvalue as it follows from the 
effective field theory to this order by inserting the above equations into eq. (|2,12|) , This is shown in 
Fig. ^for three different values of v. Clearly the approximation cannot be trusted in the region where, 
to this order, p[ (C) g° es negative. However, the comparison to the full result as it is conjectured 
from Random Matrix Theory |2U| shows that for all practical purposes even this approximation that 
keeps only the two leading terms in the expansion is sufficient when comparing with lattice data. 

From eq. 1)2.13}) and the knowledge of the two-point function eq. (|3,2j) we can also give the leading 
term in the expansion for the second Dirac eigenvalue, see Fig. [2j This 2nd eigenvalue distribution is 
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Figure 2: The distribution of the second Dirac eigenvalue for v = 0. The leading order term from eq. 
(j2.13|) (green) is plotted versus the full result (blue). 

described very accurately almost up to its maximum if we compare to Random Matrix Theory |21| , 
To get the same degree of accuracy as for the first eigenvalue we would need the three-point function, 
which has not yet been derived from the chiral Lagrangian 1 . But we hope these examples suffice to 
illustrate how individual lowest-lying eigenvalue distributions of the Dirac operator can be obtained 
from the effective Lagrangian framework. 

4 Conclusions 

We have shown that the distribution of individual Dirac operator eigenvalues can be computed solely 
from the knowledge of eigenvalue density correlations function. Since a well-established procedure ex- 
ists for deriving all /c-point spectral correlation functions from field theory, this settles the question as 
to how one can also derive from field theory individual eigenvalue distributions of the Dirac operator. 
We have illustrated the method with a case where exact analytical results already have been derived 
from field theory, namely the lowest-lying Dirac eigenvalues in QCD-like theories which undergo spon- 
taneous chiral symmetry breaking. While it is expected that all eigenvalue correlations in this regime 
are identical to those obtained from Random Matrix Theory we have not made use of this assumption. 
We have shown that even just on the basis of the knowledge, to date, of the microscopic one- and 
two-point density from effective field theory the distribution of the first eigenvalue is indistinguishable 
from the conjectured analytical result from Random Matrix Theory. The expansion can be pushed to 
any needed degree of accuracy by computing the required higher /c-point correlation functions from 
field theory. This is certainly the case for a physical number of fermions Nf. The only exception may 
be the quenched limit of Nf = 0. If this case is inherently ill-defined at the level of the expansion 
in the effective Lagrangian, not even spectral correlation functions can be reliably computed by the 
route through effective field theory. An entirely different approach could then be called for. But our 
relations between individual eigenvalue distributions and spectral correlation functions will remain 
valid. 

Acknowledgments: We wish to thank the CERN Theory Division, where part of this work was being 
done, for hospitality. The work of GA was supported by a Heisenberg fellowship of the Deutsche 
For schungsgemeinschaf t . 

1 We have checked this explicitly by using the three-point function conjectured from Random Matrix Theory. 
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A Appendix A 

Here we show that the distribution of the k-th eigenvalue eq. (|2.1U|) is normalized to unity for all k: 
dsp k {s) = (A.l) 

oo 







k \k) J dS J d ^ 1 '-- d ^ k - 1 J d\ k+ i ■ ■ .d\ N Vn{X\, . . . , Afc_i,s, Afc+i, . . . , X N ) 

k[ k )J2(-i) l [ l )L ds (i) (I ) d\ l ...d\ k - 1 d\ k+1 ...d\ N 



\J0 



x Vn(X\, ... , Afc_i, s, Afe+i, . . . , Xn) 







We have used again the identity eq. (|2,5() choosing a = J °° dX and b = dX in order to replace the 
integrals J dX. We next order the integration variables of the (N — k + 1) integrations over [s, oo), 

oo / roo\ N—k+l 

dslj ) dXk-i + i...dX N V N (Xi,...,X k -i,s,X k+ i,...,X N ) = (A.2) 

poo roo poo poo 

= {N-k + l)\ ds d\ k -i +1 dX k „i +2 ... dX N 

JO Js J A fc _; +1 JAjv_i 

x Vn(Xi, . . . , X k -i,s, Afc+i, . . . , Ajv) 

^ / f-oo\N-k+l+l 

= N _ + i + i [J J dsdX k -i + x...dXNT > N(Xi,...,X k ^ l ,s,X k+ i,...,XN), 

and then undo the ordering with the additional integration over s. Inserting eq. (|A.2|) into eq. (|A.1|) 
all the integrations over the joint probability distribution will give unity and we arrive at 

f *»W - *(^g(-l)'(*7 1 ) lf - 1E L 7TI _ , (A.3) 
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